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ABSTRACT - " 



The "perpetual salt fountain" can occur in an area of the ocean 
where the water is vertically stable. Warm, saline water overlies cold, 
relatively fresh water. A tube is placed in the water from just beneath 
the ocean surface to just beneath the halocline. Water is pumped up the 
tube to initiate the flow. As the water in the tube is warmed, it be- 
comes less dense than (and thus buoyant with respect to) the surrounding 
water, and the flow will continue after pumping is stopped. 

The effect of the vertical conduction of heat by turbulent mixing 
on the velocity of flow in the "perpetual salt fountain" is investigated. 
Temperature and salinity profiles with depth from hydrographic soundings 
are considered to determine flow velocities in the ocean. 

The calculations show that vertical heat conduction has a negligible 
effect on the rate of flow in the presence of horizontal heat conduction, 
the value of the average salinity of a column of water minus the salinity 
of the water entering the bottom of the tube is the primary factor in- 
fluencing the rate of flow, and the effect of the temperature -depth profile 
on the velocity of flow is negligible so long as initial stability re- 
quirements are met. 



u. k 

i: 



f 






TABLE OF CONTENTS 

Section Page 

1, Introduction 7 

2, Solution for a model which includes vertical heat 8 

conduction effects 

2.1 Solutions for fourth-power temperature and 10 

salinity profiles 

2.2 Solutions for linear temperature and salinity 14 

profiles 

2.3 Effects of vertical heat conduction 18 

3, Digital computer solutions 19 

4, Results 21 

4.1 Vertical heat conduction versus no vertical 21 

heat conduction 

4.2 Computer results versus Groves’ result 22 

5, Conclusions and acknowledgements 23 

Appendix 

A. Computer program for flow velocity calculations 26 

B. Computations of ^^9 and 31 

C. Wall friction calculation 32 

D. Coefficient of vertical heat conduction for turbulent 33 

flow 

E. "Fixed parameter" values 34 

Bibliography 41 



3 



LIST OF ILLUSTRATIONS 



Figure Page 

1. Initial temperature-salinity profiles constructed 35 

to approximate Groves' data 

2. Linear temperature -salinity profiles 36 

3. Actual temperature-salinity profiles 37 

4* Actual temperature-salinity profiles 38 

5. Actual temperature-salinity profiles 39 

6. Actual temperature- salinity profiles 40 

Table 

1. Results of calculations for rate of flow 17 



4 



LIST OF SYMBOLS 



w wall thickness of tube 
^ density of sea water 
c specific heat of sea water 
^ kinematic viscosity of sea water 
g acceleration due to gravity 
R Reynolds nunber 

T stress on inner surface of the tube 
K coefficient of heat conduction for turbulent flow 
v average velocity of flow in the tube 
n heat conductivity of tube material 
r inner radius of tube 
y depth below the ocean surface 
T temperature 
p pressure 

h depth of bottom of the tube 
L length of the tube 

0(y) temperature of water outside the tube at depth y 

(?(y) temperature of water inside the tube at depth y 

S(y) salinity of water outside the tube at depth y 

S(h) salinity of water at depth h, and salinity of water inside 
the tube 

# =2nL/wcrvh^ 

c*. =2n/Krw 

@ =hcvp/LK 

a=KL/cvhp 



S 



List of Symbols (Continued) 



x C$ 2 + **■)**■ 
z 

\ _ ($ 2 i- i^Y 2 - 

K j 

D = + 21 $* 

£’= 5 - 3 ] 2 i$ 3 tJ 
f- ~i(j 2C ^ + 

& = x[m] 




f \ 



1, Introduction. 



Stcmmel, Arons, and Blanchard [8] first wrote about an "oceano- 
graphic curiosity" which they called the "perpetual salt fountain." 
Consider a situation in which there is a layer of cold, relatively 
fresh water which underlies a warm, saline upper layer of water. The 
density of water is a function of temperature, salinity and pressure, 
and the density of the upper layer must be less than the density of 
the lower layer for stability. 

Now insert a tube into the water which extends from just beneath 
the surface of the water to the top of the lower layer of water. A 
power source is provided to initiate an upward flow of water. As the 

water in the lower layer is started up the tube, it is warmed by heat 

flowing into the tube from the surrounding warm water. Since the 

water in the tube is less saline than the water outside the tube, it 

becomes less dense than the outside water and is forced upward. 

The flow process requires no source of power other than that used 
to initiate the flow, and flow will continue so long as there continues 
to be a well-defined two- layer system. The principle will work for a 
flow in either direction. 

Groves [2] calculated the rate of this flew in the tube to de- 
termine what, if any, practical value this perpetual flow phenomenon 
might have. Groves’ calculation was based on a simple model where the 
water temperature was considered a constant down to the bottom of the 
tube which just extended into water of a colder temperature. Groves 
determined that the volume rate of flow for his tube would be too small 
to be of practical significance for fish farming (that is, use of this 

principle as a source for providing nutrient- rich deep water at the 
surface to increase the amount of edible fish in an area). 
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Groves' model neglected the effects of vertical heat conduction 
within the tube, which should increase the warming of the water inside 
the tube and provide a corresponding increase in rate of flow. If this 
were the case, then perhaps a larger tube could be used. 

Ihe additional heating which is a result of vertical heat con- 
duction would be sufficient to warm the increased volume of water in 
a larger tube. If the flow velocity in a tube of larger inner radius was 
approximately equal to the flow velocity in the tube of smaller inner 
radius, the tube with the larger inner radius would provide an increase 
in volume rate of flow. 

It was considered that a significant increase in the velocity of 
the water in the tube might increase the volune of flow sufficiently to 
warrant further investigation of this process and its practical appli- 
cations. Hie effects of wall friction on the temperature of the water 
inside the tube were considered, but the result was so small that the 
effect was considered negligible. Calculations of temperature change 
due to this stress are contained in Appendix C. 

2. Solution for a model which includes the effects of vertical 
heat conduction . 

In order to see the effect of vertical heat conduction on the 
rate of flow of the water inside the tube, let us consider the same 
model that Groves worked with, but include vertical heat conduction. 
Assume that the water density is a linear function of temperature 
and salinity, and that the water is incompressible. Therefore, the rate 
of change of density with temperature, , and the rate of change 

of density with salinity, , are constant, a valid assumption 

within the small range of temperature and salinity considered. Then, 

? = + & (*4e\ + S (*%s) a 0) 
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where 9 is the water density and 0 and S respectively are the temp- 
erature and salinity of the water. The heat flow into the tube is 



calculated by the equation 

S (s** 1 d if) = /« 2 ( -^K 



w 



Kir ifz 



( 2 ) 



where D/Dt designates the material derivative following the mean 
motion of the water, q is the heat per unit mass, n is the heat con- 
ductivity of tube material ( say, copper ) , r is the inner radius of 
the tube, 0 (y) is the temperature of the water in the tube, w is the 
wall thickness of the tube, and K is the coefficient of vertical heat 
conduction for turbulent flow calculated in Appendix D. At the sea 
surface y=0, and y is positive downward. If we assume a steady state, 

00 „ c^P. jjj _ _ <^<P 

Dt " " -Vr (3) 

since the local change ( The average flow velocity is v.) 



(4) 



Combining the constant terms, equation (2) becomes 

+ fc 1$ - c* Cp = - oC© 

where (^hcv^/LK and <X=2n/Krw. The specific heat of sea water is c, h is 
the depth of the bottom of the tube, and L is the length of the tube. 

The general solution of the homogeneous equation associated 
with (4) is 



where Jj and j£ are undetermined constant coefficients, 

v . + O* 

. ~ 2 
and , 

Az 2 • 



C5) 

( 6 ) 

(7) 
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2.1. Solutions for fourth-power temperature and salinity profiles . 

In order to find a particular solution to equation (4) , a fourth- 
power temperature profile was assumed: 

- A' . ( 8 ) 

This function was chosen because it can be easily fit to actual 
temperature -depth curves. The value of sea-surface temperature can 
be inserted for A, and the value (T s - T^J/h^ inserted for C. The 
complete solution for equation (4) is _ 

+ ^e - * 2 * + A - -^-\2 l IC‘< 2 + 72tfc*+2‘1p' , Cj 

' is f Mpc * + 21 *] if - 52 [l2C* + 1 2 fr] 

- dr foe} I 3 - C a y . 

The water entering the bottom of the tube is at the same temper- 
ature as the surrounding water at that depth; therefore, 0(h) s 9(h). 

For steady flow, the water at the top of the tube is assumed to be at 
the same temperature as the water surrounding it. This is a reasonable 
assumption since we are considering vertical heat conduction, and the 
water at the top of the tube undergoes vertical turbulent mixing with 
the surrounding water. Thus, the other boundary condition is 
0(0) =9(0), Only upward flow is considered in this paper. 

Applying the boundary conditions to equation (9) and solving for 



and ^ p 

<P(i) =• 



p(i-e. A * h )+ Eh + Fh 2 



A,h 



- t' 2 

~P(V‘ K -Q - Eh - Fh* e A 2 f 



(10) 



+ A- D - fiY * f - G y 3 - 0 

The constants D, E, F, and G are functions ofo<, @ , and C for a given 
velocity: 



10 



d = 3 [zic^ 2 1- 72? J c* + 

E = 53 j VS’ptot + Zi^cji 

p . ^[ 1 * c=t ■*• 12 e 2 c] 5 
<S- 5rCv?c] . 



Now assume that the water is removed from the top of the tube as 
fast as it flows up so that no additional pressure head is formed. 

Then the pressure difference driving the flow is the hydrostatic pres- 
sure outside the tube at y=h, less the hydrostatic pressure in the tube 
at this depth. 



(ID 



( 12 ) 



Equation (1) gives an internal pressure difference of 

A f^ 1 

Jo 

and an external pressure difference of 

foM- fo Co) = fjj?. + (x4 s(n) + (H) o ©6)] d y # 

Since both pressures are atmospheric at the surface, subtraction of 
equation (12) from equation (11) gives the pressure which drives the 
flow: h 

A-p = + j©6D‘^)jieJ dy . 

Dfensity is constant to within one per cent, and from the continuity 
equation, ^v=constant. This implies that viconstant and, in fact, 
will be assumed constant. Then the wall friction, T’, is also constant: 



(13) 



r= g f(R) ?\r 2 . 

Here, f(R) is called the "resistance coefficient," (Rouse [5]). Due 
to the small value of viscosity for sea water, the flow will be as- 
sumed to be turbulent, and the relationship 
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(15) 



■F(r) = . 3H,i/R /lf 



(16) 



for smooth tubes will be used. Here, R is the Reynolds Number, 

ft = 2\T/T / 

where is the kinematic viscosity of sea water. 

Since the water in the tube is not accelerating, the driving 
force due to pressure must equal the retarding force due to wall 
friction: 

/\ 2 LT (17) 

- — * 



Combining equations (14) and (17) gives 



L f(R) 



U 



~ A -p - o . 



(18) 



Solving explicitly forZl-pin equation (13) necessitates prescribing 
a salinity- depth curve. A fourth-power salinity profile 

S(^) = Q - a 9 ) 

gives the salinity curve with depth the same physically reasonable 
form as the temperature curve. The surface salinity is Q and 

O A n.4 



( 20 ) 



Equation (13) becomes 


l (l-e A * h )- 


u-4 


L S J d L ^ 1^2 J 



*- - 1 U - i 

_ \(eu+ - A 2 (i- 

~ he L J * 



Equation (20) is difficult to evaluate in its present form, but 
there are several approximations which can be made to simplify’ the ex- 
pression without sacrificing accuracy. 
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First, assune that the rate of flow i^ill be less than ten meters 
per second. Calculated values were mudi less than this value, and a 
velocity of ten meters per second would require an extremely large 
pressure force to drive the flow. 

Second, assume that the tube thickness is very much less than the 
inner radius of the tube, A thin tube is required to allow the water 
in the tube to reach approximately the same temperature as the water 
outside the tube as it flows upward. 

Third, assune terms that are 10"^ (order of magnitude) as large 
as the dominant teims are negligible. 

From equation (5) , Aj , is a function of «<• and and will always 
be a negative number. If h is greater than (say) 50 meters, (it will 
usually be much greater than this), A|h is a large negative number, 

X h 

and e A{ is negligible. 

The constants D, E, F, and G are such small numbers that, even 

2 3 

when multiplied by 1, h, h , and h respectively, they are still negli- 
gible compared to other terms. 

Equation (20) becomes 

Av = ^Tsfe-ir- 5 ^] (2D 

and combining equations (18) and (21) gives 

£$]- ,hjf [«- T - ■&>]■ ™ 

This is a closed foim of the solution for the flow velocity which in- 
cludes the effects of vertical heat conduction. 

If vertical heat conduction is not considered, the governing 
differential equation (2) simplifies to 
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(23) 




= m 



2TTsrdy 



(e-d)) 



w 



which has the solution 




where & =2nL/wcvrh^> . 

Solving equation (24) for the fourth-power temperature profile gives 



p/Ck 3 , /2Ck 2 . l)~] 

" £ L~*~ J ' 

Using the fourth-power salinity profile and equation (25) for 0(y), 
equation (13) becomes 



(25) 




(26) 



where negligibly small terms have been ignored. The solution for fourth- 
power temperature and salinity profiles when vertical heat conduction is 
not considered is 




PK 







(27) 



Equations (22) and (27) were solved using the temperature and sali- 
nity profiles pictured in figure 1, using identical tube dimensions. 

The flow velocity which was influenced by vertical heat conduction was less 
than one- tenth of one per cent greater than the flow velocity calculated 
ignoring vertical heat conduction (Table 1,), 

2.2 Solutions for linear temperature and salinity profiles . 

In order to consider the effect of different temperature and sali- 
nity profiles, linear profiles were assumed: 
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and 



0(%) = A\- Ny, 
S(y) - ~ ^ * 



( 28 ) 



(29) 



The solution to equation (4) then becomes 




where X=N^. 



Solving equation (13) using the temperature profile of equation (30) 



gives 



A-p =?bf s 



Q- 




6W<) 

A i 



(31) 



where negligibly small terms have been ignored. The solution for linear 
temperature and salinity profiles, including vertical heat conduction 




Assuming linear temperature and salinity profiles, equation (24) becomes 



= ] + fA-Ahf (53) 

and the solution to equation (13) is 

where negligibly small terms have been ignored. Thus, the solution for 
flow velocity for linear temperature and salinity profiles where verti- 
cal heat conduction is ignored is 
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V 




* (35) 



Equations (32) and (35) were solved using the temperature and sali- 
nity profiles pictured in Figure 2, with identical tube dimensions, Again, 
the rate of flow influenced by vertical heat conduction was less than one- 
tenth of one percent greater than the rate of flow calculated ignoring 
vertical heat conduction, (Table 1,), 

It is interesting to compare the actual calculations of flow veloci- 
ty for the profiles in Figure 1 with those related to Figure 2, Vertical 
heat conduction was included, but this is not pertinent to the discussion 
below. The velocity solutions, equations (22), (25), (32) and (34) all 



between the surface and the depth y=h, S, less the salinity at y=h, S(h), 
In making these calculations it was noted that the main factor causing 
the one centimeter per second difference in rate of flow between the pro- 
files in Figures 1 and 2 was the value of the S-S(h) term. 

Calculating the value of S-S(h) for the salinity profiles in Figure 
1 and in Figure 2 showed that the two values were almost identical. 
However, the value of IT obtained by applying the fourth-power salinity 
profile approximation to the salinity curve of Figure 1 gave an erroneous 
value for S-S(h), The error in the calculation of S by substituting an 
approximate function for the actual salinity curve caused a marked error 
in flow velocity result. This error was later verified by the computer 
program determination of flow velocity which required no approximation of 
"s’. This indicates the great importance of the salinity profile which 
will be discussed in more detail below. 



contain the term 




which is the average salinity of the water 
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The flow rates for the 9-S profiles contained in Figures 1 and 2 were 
also calculated utilizing a digital computer (see below) , Here the ef- 
fects of vertical heat conduction were not included. The rate of flow in 
each case was within one per cent of the other, and was essentially the 
same as the flow rate calculated by hand from equation (32) , 

TABLE 1 

RESULTS OF CALCULATIONS FOR RATE OF FLOW 



Case 


0(y) $ S(y) 


"Fixed # 


Hand ## 


Computer ## 


No. 


Profiles 


Parameters" 


Calculations 


Calculations 


1 


Figure 1 


Approximate 




17.32 


2 


Figure 2 


Approximate 


17.36 


17.30 


3 


Figure 1 


Approximate 

** 


16.32* 




4 


Figure 2 


Approximate 


28.52 


27.88 


5 


Figure 1 


Precise 




17.40 


6 


Figure 3 


Precise 




22.68 


7 


Figure 4 


Precise 




26.60 


8 


Figure 5 


Precise 




20.89 


9 


Figure 6 


Precise 




8.03 


10 


Figure 5 


Precise*** 




30.94 


11 


Figure 3 


Precise*** 




33.59 



* The flow velocity result is for the hand calculation of the fourth- 
power approximation of temperature and salinity profiles illustrated in 
Figure 1, 

** The inner radius of the tube was increased to 20 centimeters and the 
thickness of the tube increased to two centimeters. 

*** The tube was inserted vertically so that the length of the tube equals 
the depth of the bottan of the tube. 

# Appendix E contains "fixed parameter" values. Except where otherwise 
noted, the inner radius of the tube is 0.1 meters and the wall thickness 
is 0.002 meters. 

## Flow velocity results are in centimeters per second. 
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Marked differences in rate of flow which occur when the tube is 
inserted into water which has different temperature and salinity pro- 
files are a direct result of the difference in the salinity profile from 
one location to another. The temperature profiles did not significantly 
influence the rate of flow within the small but realistic range of 
temperature and salinity used by the author. 

2,3 Effects of vertical heat conduction . 

The importance of vertical heat conduction was investigated in 
another manner. The wall thickness of the tube was increased by a factor 
of ten. The flow velocity decreased, but the decrease was less than one 
per cent. Increasing the inner radius of the tube from ten centimeters 
to 20 centimeters and retaining the thicker-walled tube resulted in a 
significant increase in rate of flow of 371. 

Comparison of the results of the above calculations indicated that 
the vertical heat conduction due to vertical turbulent mixing was suf- 
ficient to warm the water in the tube. In order to determine the extent 
of the wanning due to vertical heat conduction, the velocity of flow was 
then determined for the situation in which no horizontal heat conduction 
was allowed. In this case, the temperature in the tube is a function 
only of vertical heat conduction. The governing differential equation (2) 
simplifies to 




(36) 



which has the solution 




where a=KL/cvh^> , 
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Solving equation (13) using equation (26) and linear temperature and 
salinity profiles , 

where negligibly small terms have been omitted. 

The resulting flow velocity is nearly identical to the flow velocity 
calculated using the same tube, but allowing horizontal heat conduction. 

The process of heat conduction in the vertical by turbulent mixing would 
require a longer period of initial pumping to reach a steady state con- 
dition, but the results indicate that the vertical heat conduction pro- 
cess is sufficient to heat the water inside the tube. 

The initial tube placement was at an angle of depression of 30 degrees 
to the water surface. This was done so that the water in the tube would 
require a longer time interval to reach the surface and undergo increased 
warming. Since the water in the tube is sufficiently warmed by vertical 
heat conduction to be less dense than the water surrounding it outside the 
tube, the tube placement was changed to a vertical position. Hie decrease 
in tube length resulted in a corresponding decrease in the effect of wall 
friction and a 481 increase in rate of flow. Such an increase is certainly 
significant and will increase the volume rate of flow. 

3. Computer solutions. 

In order to determine the rate of flow for a "perpetual salt fountain" 
situated in the ocean, actual oceanographic temperature and salinity data 
for a given area must be known. A program to obtain flow velocity so- 
lutions which neglected vertical heat conduction was written for a CDC 
1604 digital computer using FORTRAN 60 programming. Using this program, 
flow velocity was easily determined for any area in which temperature and 
salinity data were available. 






\f ‘ijb.loStlD 
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Temperature and salinity values for each one-meter interval over the 
entire length of the tube, tube dimensions , and "fixed parameters" are en- 
tered on data cards in the program, "Fixed parameters" include kinematic 
viscosity of sea water (^ ) , density of sea water (^>), specific heat of 
sea water (c) , acceleration of gravity (g) , rate of change in density with 
change in salinity C^/^S) and rate of change in density with change in 
temperature » Appendix E contains the values used for "fixed 

parameters," 

Subroutine INTE computes the temperature profile in the tube, and 
subroutine DELTP computes the pressure force which drives the flow. Sub- 
routine INTE provided some interesting problems due to the size of the 
exponents of the exponential functions, A discussion of these problems, 
and the solution to them is included with the computer program in 
Appendix A, Subroutine DELTP is a trapezoidal integration which is as 
accurate as can be obtained from the small number of data points gene- 
rally given in a hydrographic sounding. 

Initial data used in the program was set up to treat the same case 

as was calculated by Groves, Temperature and salinity profiles were 
devised so that values of 9(0) -0(h), "S-S(h), and the "fixed parameters" 
would be identical to those used by Groves, Additional test cases were 

run to compare with the results of hand calculations. 

Values for the constant teims^eand J %s were then calculated from 

the equation of state given by Eckart [1] (see Appendix B) and more pre- 
cise values for the other "fixed parameters" were obtained and used with 
data from actual hydrographic soundings. Several cases of actual data 
were employed to get an idea of the range of values of flow velocity which 
could be expected in the ocean. The results are contained in Table 1, and 
the temperature-depth and salinity- depth curves are pictured in Figures 1-6. 
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In the actual ocean data cases, the depth to which the tube was ex- 
tended was determined by the depth of the halocline. Hie tube was always 
extended into the water just below a negative halocline, 

4. Results , 

4.1 Vertical heat conduction versus no vertical heat conduction . 

Initial flow velocity calculations were made by hand using the initial 
values of the "fixed parameters" contained in Appendix E. The rate of 
flow for identical tube dimensions was calculated for 6-S profiles con- 
tained in Figures 1 and 2. 

The calculation of flow rate was made twice for the linear 9-S pro- 
files of Figure 2, once including the effects of vertical heat conduction 
and once neglecting the effects of vertical heat conduction. In each 
case, inclusion of vertical heat conduction had a negligible effect on 
the rate of flow. 

Flow velocity calculations were then made for the 9-S profiles of 
Figure 1 using the fourth-power temperature and salinity profiles. These 
functions were selected because they approximated the physical appearance 
of the curves in Figure 1 and because they were easily integrable. Flow 
rate solutions for the vertical heat conduction and no vertical heat con- 
duction equations were again identical to four significant figures. 

Flow velocity was then calculated for a situation in which no hori- 
zontal conduction of heat was allowed, so that only vertical heat con- 
duction was possible. Hie rate of flow for this calculation was nearly 
identical to the rate of flow calculated for two other cases which used 
the same tube dimensions, "fixed parameter" values, and 9-S profiles. 

Hie first allowed both horizontal and vertical heat conduction, and the 
second allowed cnly horizontal heat conduction. 



21 



On the basis of these calculations, it was concluded that the effect 
of vertical heat conduction on the rate of flow of the water inside the 
tube was negligible in the presence of horizontal conduction. Thus it 
was decided to make use of the computer program which did not include 
the effects of vertical heat conduction for the remaining calculations, 

4,2 Computer results versus Groves* result . 

Utilizing the computer program allowed exact values of temperature 
and salinity at each depth to be used. Temperature and salinity values 
from Figures 1 and 2 were used in the program with a resulting rate of 
flow of 17.30 centimeters per second (cps) for Figure 1 profiles and 
17,32 cps for Figure 2 profiles. The 17,32 cps flow velocity associated 
with the profiles of Figure 2 was in close agreement with the rate of 
flow calculated by hand for the same profiles. The 17.30 cps rate of 
flow calculated for the profiles in Figure 1 was much greater than the 
rate of flow calculated by hand for the same profiles. Investigation of 
the hand calculation revealed that the"S-S(h) term calculated by hand 
differed significantly from the actual value of S-S(h) obtained from 
Figure 1. This difference was attributed to use of the fourth-power 
salinity profile as an approximation of the actual salinity profile 
when making hand calculations. S-S(h) calculated from the curve of the 
fourth-power salinity profile is significantly different from the true 
value of S-S(h) in Figure 1. The influence of the approximation to the 
form of the temperature profile, however, was determined to be negligible. 

The profiles of temperature and salinity in Figures 1 and 2 were 
constructed so that the values of S-S(h) and 9(0) -0(h) would be identical 
to those selected by Groves, Tube dimensions and "fixed parameter" 
values were also those used by Groves. Groves' calculated rate of flow 
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was 17,50 cps as compared with a rate of flow of 17,32 cps obtained from 
the computer program. The difference in flow velocities can be attribu- 
ted to a rounding-off procedure used by Groves to simplify calculations. 
The values of and S calculated in Appendix B were then 

inserted into the original program along with more precise values of the 
other "fixed parameters" to determine what effects the more precise values 
would have on the rate of flow. The resulting velocity increase of eight 
centimeters per second was not significantly greater, but the more pre- 
cise values for the 'fixed parameters" were retained for the remaining 
temperature and salinity profiles. 

Temperature and salinity data from actual hydrographic soundings were 
also used. Figures 3 through 6 illustrate the actual profiles used. It 
was verified that areas of most pronounced negative haloclines produced the 
greatest rates of flow, 

A marked increase in rate of flow was noted when the tube was in- 
serted vertically instead of at a depression angle of 30 degrees to the 
ocean surface. The 481 increase in speed of flow will provide an equal 
increase in volume rate of flow. 

Rate of flow was also increased by increasing the inner radius of the 
tube. This is doubly significant since an increase in flow velocity in- 
creases volume rate of flow, and increased tube radius will also increase 
tlie volume rate of flow, 

5. Conclusions and acknowledgements . 

The test cases discussed in this paper have demonstrated some interest 
ing features. First, the temperature -depth structure is of relatively 
little importance. The surface layer must be sufficiently warm to in- 
sure stability, but the magnitude of the temperature change between the 
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surface and the bottom of the tube has a negligible effect on the velo- 
city of flow. 

Second, the primary parameter influencing the rate of flow of water 
in the tube is the difference between the average salinity in the layer 
through which the tube extends, and the salinity of the water entering 
the bottom of the tube,, To calculate the rate of flow for a given case 
to within one per cent of its precise value, a salinity-depth curve must 
be known, A simple area-averaging process can be used to determine S, 
the average salinity, and S(h) is known. Then, 

/fasao'* gQf] ±lh%[s-S(h)j. (s9) 

Third, equation (39) gives excellent results because the effects 
of vertical heat conduction heat the water in the tube sufficiently to 
make it buoyant. Although vertical heat conduction does not directly 
increase the velocity of flow, the warming effects allow the substi- 
tution of a poorer heat-conducting material for the copper tube. A 
metal which is resistant to sea water corrosion would be a better material 
for the tube. In addition, a relatively thick tube with a larger inner 
radius may be used. This will increase the volume rate of flow and 
provide a more durable tube. 

For example, using data from the equatorial Atlantic Ocean region 
in case #11 (Figure 3), the volume rate of flow is 10.5 liters per second. 
This is a 90% increase in volume rate of flow over that calculated by 
Groves. When the inner radius of the tube was enlarged in case #4 
(Figure 2), the volume rate of flow increased to 34,6 liters per second. 
This is an increase of 530%, which is definitely significant, 

a 

It is therefore concluded that the "perpetual salt fountain" is more 
than an "oceanographic curiosity," Upward flow of nutrient-rich water 
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from well below the surface in sufficient vclune could make a fertile 
fishing area out of an area which is, at present, barren. A relatively 
inexpensive and simple method for providing this nutrient-rich water, 
such as the method described in this paper, is worthy of additional 
experimentation . 

The author gratefully acknowledges the guidance, encouraging 
assistance and critical comments of Assistant Professor T. Green. My 
thanks also go to Mr. R. D. Brunell of the United States Naval Post- 
graduate School’s computer facility staff for his valuable programming 
assistance. 
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APPENDIX A 



. / 

COMPUTER PROGRAM FOR FLOW VELOCITY CALCULATIONS 



. . JOBO 532 FM,H INMAN *K.G. »FLOWUP 
PROGRAM FLOWUP 

DIMENSION T ( AO 1 ) * S(401), PHI(401) 

COMMON NVT »H*GAMMA»DELH»PHI » S » T »SALTH»THETH»G»DELTP»CHDS»CHDT 
READ 50 » NVT 
50 FORMAT (1215) 

READ 10* ( T ( I ) » S ( I ) » I = 1 »NVT ) 

10 FORMAT (6E12.6) 

READ 10* R*E»CON*PL*W.D.SH*DEPTH.VO»EP1*THETH 
READ 10* H»G,DELH,SALTH,CHDS»CHDT 
PRINT 12* R»E»CON»PL»W»D 
PRINT 12* SH*DEPTH»V0»EP1*THETH 
PRINT 12, H»G»DELH»SALTH*CHDS»CHDT 
12 FORMAT (6E12.6) 

V = V 

19 V=V+0.01 
RE=2.0*V*R/E 

GAMMA=4 • 0*CON*PL/ ( W*D*SH*E*RE*DEPTH ) 

CALL INTE 
CALL FINAL 

I F ( DEL TP-0 • 0 ) 55*55*56 

55 V=V-0.01 
GO TO 57 

56 TAU=DELTP*R/( 2.0*PL) 

I F ( RE-2000. 0 ) 20,20*21 

20 FUNC=64.0/RE 
GO TO 22 

21 FUNC=0.3164/( (RE)**( 1.0/4. 0) ) 

22 VEL=(8.0*TAU/(FUNC*D) )**( 1.0/2. 0) I 

11 FORMAT ( 3 OH THE UPWARD FLOW VELOCITY IS 2E20.6) 

IF(VEL-V) 55*24,19 

57 V=V+0.001 
RE=2.0*V*R/E 

GAMMA=4«0*CON*PL/(W*D*SH*E*RE*DEPTH) 

CALL INTE 
CALL FINAL 

I F ( DELTP-0 • 0 ) 40*40,41 

40 V=V-0.001 
GO TO 42 

41 TAU=DELTP*R/(2.0*PL) 

I F ( RE-2000 • 0 ) 43,43,44 

43 FUNC=64.0/RE 
GO TO 45 

44 FUNC = 0.3164/ ( (RE)**(0.25) ) 

45 VEL=(8.0*TAU/(FUNC*D) )**(0.5) 

IF(VEL-V) 40*24,57 

42 V=V+0.0001 
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RE=2.0*V*R/E i 

GAMMA=4.0*CON*PL/( W*D*SH*E*RE*DEPTH) 

CALL INTE 
CALL FINAL 

I F ( DELTP-0 • 0 ) 46*46*47 

46 V=V-0.0001 
GO TO 70 

47 TAU=DELTP*R/(2.0*PL) 

I F ( RE-2Q00 • 0 ) 71,71,72 

71 FUNC=64.0/RE 
GO TO 73 

72 FUNC=0. 3164/ ( (RE) **( 1*0/4. 0 ) .) , 

73 VEL=(8.0*TAU/(FIJNC*D) )**( 1. 042.0) 

IF(VEL-V) 46,24,42 

70 V=V+0. 00001 v ’ 

RE=2.0*V*R/E 
PRINT 60, RE 

60 FORMAT (25H THE REYNOLDS- NUMBER IS E20.6) 
GAMMA=4»0*CON*PL/(W*D*SH*E*RE*DEPTH) 

PRINT 61, GAMMA 

61 FORMAT ( 14H GAMMA EQUALS E20.6) 

CALL INTE 

CALL FINAL 
PRINT 63, DELTP 

63 FORMAT ( 14H DELTP EQUALS E20.6) 

TAU=DELTP*R/ ( 2.0*PL ) 

PRINT 62, TAU 

62 FORMAT (12H TAU EQUALS E20.6) 

I F ( RE-2000, 0 ) 74,74,75 

74 FUNC = 64 »0/RE 
GO TO 76 

75 FUNC=0.3164/ ( (RE)**( 1.0/4. 0) ) 

76 VEL=(8.0*TAU/ (FUNC*D) )**( 1.0/2. 0) 

PRINT 11, VEL,V 

DIFF=ABSF( ( VEL-V) /V) 

IF(DIFF-EPl) 24,24,70 
24 STOP 
END 

SUBROUTINE INTE 

DIMENSION T ( 401 ) , S(401)» PHI(401) 

COMMON NVT, H, GAMMA, DELH, PHI ,S,T, SALTH,THETH,G, DELTP, CHDS,CHDT 
DO 10 J=1 ,NVT 
PHI ( J)=0.0 
10 CONTINUE 
I J=NVT 

PHI ( IJ)=THETH 
BIG=GAMMA+1.0 
BIGEX=EXPF( GAMMA ) 

BNE=GAMMA-1.0 
20 I J= I J- 1 
K= I J+l 

X I = ( T ( IJ)-(T(K)*BIG) >/ ( GAMMA*B I GEX ) 
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XJ=( (T( IJ)*8NE)+T(K) ) /GAMMA 
XK=PHI ( K ) /8IGEX 
PHI ( I J ) =,X I + X J + X K. 

IF(.IJ-l) 30,30*20 
30 RETURN 
END 

SUBROUTINE FINAL 

DIMENSION T ( 40 1 ) * S(401), PHI(401) 

COMMON NVT , H *GAMMA , DELH , PH I *S , T , SALTH . THETH. G ♦ DEL TP ♦ CHDS ♦ CHDT 
TOT AL = 0 «0 ( 

1 = 1 

PRESS= ( S ( I ) -SALTH ) *CHDS+( T ( I ) -PHI ( I ) )*CHDT 
40 1=1+1 

SUMP= ( ( S ( I ) -SALTH ) *CHDS ) + ( ( T < I )-PHI ( I ) ) *CHDT ) 

TOM=DELH*( PRESS+0. 5* (SUMP-PRESS) ) 

PRESS=SUMP 
TOTAL=TOTAL+TOM 
IF(NVT-I) 50*50*40 
50 DELTP = G*TOT AL 
RETURN 
END 
END 






t 
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The computer program is an iteration process for obtaining the 
velocity of the water in the tube. It was written to solve actual 



oceanographic cases using the form of the solution given by Groves [2] . 
Subroutine DELTP is utilized to compute the pressure integral, 



The subroutine is a trapezoidal integration which assumes a linear 
relationship between successive points. No greater accuracy is desirable 
because hydrographic soundings contain salinity and temperature data for 
relatively few depths, and the linear relationship is generally assuned 
when plotting the original data. 

Subroutine INTE is designed to compute the temperature of the water 
at each one-meter interval of depth inside the tube using 



Equation (A2) relates water temperature in the tube to depth, but it 
cannot be solved in this form on the computer since the capacity of the 
CDC 1604 is exceeded by exponential functions with exponents of magni- 
tude greater than 709. A different form of equation (A2) was necessary 
for the computer program. ^ 





(A2) 



Let 



and 




Then 




(AS) 



Jn 4 - 7 V» l.m 1 4 - 




(A6) 
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and then from equations (A3) and (A6) 

XG-aK) + e' tAh i(i). <"> 



We now assume that 9(T) , where T is a "dummy" variable of integration,, 
varies linearly between consecutive one-meter intervals, and limit the 
interval, Ah, to one meter. 



Combining equations (A5) , (A7) , and (A8) 

<p6j-av>) = -j^-er'tfeC^Ci-n ) ' e (i|-AW)j 

+ [a(^) + £ V6). tA9) 

Since 0(h)=9(h) is a boundary condition, the first 0(y-Ah)will be cal- 
culated for a depth of (h-1) meters. This integration subroutine will 
function so long as &< 709. 
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APPENDIX B 

COMPUTATIONS OF ^9 AND ^%S 



The values for the rate of change of density with changing tempera- 
ture and salinity were assumed constant in the solution of the pressure 
integral. Additional investigation of the value of these constants was 
considered important. 

Assuming that the effects of pressure are negligible, 

. < B1 > 

Following Eckart [1 ] t 

(v + Vo)(A 0 -A ) = K w 

where 

A = /O* (l- y ) (B3) 

and 



-po = S210 +• 3PT -h O. 3 7<TT /2 +■ 3 5 

\( = jj i. ws * a nzsT-o.oooi is r x -(d.o3%o+o. 

Z \ 0 = constcin'i: = 3020 , 



(B4) 

oooir)sJ x/0 \ K) 

m 



Combining equations (B2) and (B3) gives 




Taking the partial derivatives of ^ with respect to 0 and S in 
equation (B7) , and evaluating the result at an average temperature, 
0=15°c, an average salinity, S=34°/oo, and a zero pressure, gives 
^?/6s =0.77667kg/m 3 *o/ooand ‘tyg =0. 19987kg/m 3,0 C. 
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APPENDIX C 



WALL FRICTION CALCULATION 



It was considered that the effects of the stress on the water by 
the wall of the tube might cause significant heating of the water. The 
work done on the water per unit mass per time is 



T 2TT /f dy sT 


,2 TV- 


( joules A 




^ IT 




(^kg°sec y. 


(Cl) 


The work per unit mass is 

ZT V _ Lj 

(p 'f \r 


2Tb 

9* 


• 

in 

<1> 

r-H bfl 

8* 

•rn 


(C2) 



Dividing by the specific heat, c, and converting joules to calories using 
0,2389 calories/ joule at 15°C, 



-^■(o.23P?) = AT (°C). CC3) 

Applying the data used in the solution of case #1, the temperature change 
of the water in the tube due to wall friction is 0,00025°C for a velocity 
of 17,32 cps and a pipe length of 600 meters. This change in temperature 
is negligible. 
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APPENDIX D 



COEFFICIENT OF VERTICAL HEAT 
CONDUCTION FOR TURBULENT FLOW 



The value for the coefficient of vertical heat conduction for turbulent 
flow, K, was investigated. From Schlichting, [6], the coefficient of eddy 
conductivity of heat, A^, is at least as large as the coefficient of eddy 
viscosity, A r , over the entire pipe. We assume the turbulence to be ap- 
proximately isotropic, so that turbulent transfer along the pipe is the 
same order of magnitude as that in any cross section. 

A good estimate of the effective value of Aq can then be obtained by 
calculating an average of A-j-over a cross section. This is easily done 
using mixing length theory. The result is 



v/T 






= A r = k 2 M*2 = 9 ^ u *_£. 



7I> 2 



(Dl) 



where r is the interior radius of the tube, z is measured inward from the 
wall, u* is the "shear velocity" (see Schlichting), and k=0,4. 

Then, 



K - Af c p ♦ 



(D2) 



This has a typical value of 5.0 cal./m* °Osec, for the rates of flow 
and tube size considered. 
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APPENDIX E 



"FIXED PARAMETER" VALUES 

The expression "fixed parameters" refers to values for kinematic 
viscosity of sea water (*4), conductivity of tube material ( n) , densi- 
ty of sea water (£>), specific heat of sea water ( c ), acceleration of 
gravity ( g ), rate of change of density with changing salinity (^P/^ S) 
and rate of change of density with changing temperature (^P/£)©) , Two 
sets of values for the "fixed parameters" have been used. One set con- 
tains rounded- off values for hand calculations, and the other set contains 
more precise values for the computer program. One could then determine the 



effect on rate 


of flow of changes 


in the values of the 


"fixed parameters," 




APPROXIMATE 


PRECISE 


UNITS 




10‘ 6 


1.22X10" 6 


m 2 /sec 


m (copper) 


10 2 


io 2 


cal/m*°c,sec 


9 


10 3 


1.026X10 3 


kgm/m 3 


C 


10 3 


9.33X10 2 


cal/kgm«°C 


* 


10 


9.8 


m/sec 2 


V<)S 


0.7 


0.77667 


kgm/m 3 




-0.2 


-0.19987 


kgm/m 3 *°C 
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INITIAL TEMPERATURE-SALINITY PROFILES 
CONSTRUCTED TO APPROXIMATE GROVES 1 DA.lTA 




Figure 1. 
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LINEAR TEMPERATURE-SALINITY PROFILES 
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Figure 2 



ACTUAL TEMPERATURE- SALINITY PROFILES 




Figure 3. 



37 



ACTUAL TEMPERATURE-SALINITY PROFILES 




Figure 4. 



38 



ACTUAL TEMPERATURE-SALINITY PROFILES 
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Figure 5 . 
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ACTUAL TEMPERATURE-SALINITY PROFILES 
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